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Abstract. We have developed a method for selecting the most 
luminous galaxies detected by IRAS based on their extreme 
values of R, the ratio of 60 /im and £>-band luminosity. These 
objects have optical counterparts that are close to or below 
the limits of Schmidt surveys. We have tested our method on 
a 1079 deg 2 region of sky, where we have selected a sample 
of IRAS sources with 60 /im flux densities greater than 0.2 Jy, 
corresponding to a redshift limit z ~ 1 for objects with far- 
IR luminosities of 10 13 L Q . Optical identifications for these 
were obtained from the UK Schmidt Telescope plates, using 
the likelihood ratio method. Optical spectroscopy has been car- 
ried out to reliably identify and measure the redshifts of six 
objects with very faint optical counterparts, which are the only 
objects with R > 100 in the sample. One object is a hyper- 
luminous infrared galaxy (HyLIG) at z = 0.834. Of the re- 
maining, fainter objects, five are ultraluminous infrared galax- 
ies (ULIGs) with a mean redshift of 0.45, higher than the high- 
est known redshift of any non-hyperluminous ULIG prior to 
this study. High excitation lines reveal the presence of an ac- 
tive nucleus in the HyLIG, just as in the other known infrared- 
selected HyLIGs. In contrast, no high excitation lines are found 
in the non-hyperluminous ULIGs. We discuss the implications 
of our results for the number density of HyLIGs at z < 1 and 
for the evolution of the infrared galaxy population out to this 
redshift, and show that substantial evolution is indicated. Our 
selection method is robust against the presence of gravitational 
lensing if the optical and infrared magnification factors are sim- 
ilar, and we suggest a way of using it to select candidate gravi- 
tationally lensed infrared galaxies. 

Key words: galaxies: distances and redshifts - galaxies: evo- 
lution - galaxies: starburst - infrared: galaxies 



1. Introduction 

One of the most fundamental results of the IRAS survey was 
the discovery of a class of galaxies that emit the bulk of their 
energy at far-infrared wavelengths and that have 8 — 1000 /im 
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luminosities > 10 L Q (for H = 75kms Mpc and 
go = 0.1, as assumed throughout this paper). As shown by 
Soifer et al. (1986, 1987), these ultraluminous infrared galax- 
ies (ULIGs) are a very significant population in the local uni- 
verse, dominating the high luminosity end of the local luminos- 
ity function, and outnumbering local optically selected quasars 
by a factor of at least 2 (see Sanders & Mirabel 1996 and refer- 
ences therein). More recently, a number of objects with Lir > 
10 13 Lq have been identified, for which the term hyperlumi- 
nous infrared galaxies (HyLIGs) has been proposed. The first 
object of this type to be discovered in an infrared-selected sam- 
ple was IRAS F10214+4724 at z = 2.28 (Rowan-Robinson et 
al. 1991a). Although this object is now known to be gravita- 
tionally lensed (Broadhurst & Lehar 1995; Close et al. 1995; 
Serjeant et al. 1995; Eisenhardt et al. 1996), it is a HyLIG 
even after the gravitational magnification has been accounted 
for (Downes et al. 1995). Furthermore, of the now more than 
30 known objects in this class, most do not show any evidence 
of gravitational lensing. 

The nature of the energy source powering the high bolomet- 
ric luminosities of ULIGs and HyLIGs is a subject of intense 
debate (see e.g., Rowan-Robinson 1996 for a recent review). 
Since ULIGs and HyLIGs have luminosities in the range of 
the most powerful active galactic nuclei (AGNs), it has been 
proposed that they are powered by dust-embedded AGNs. Al- 
ternatively, the far-infrared (FIR) luminosity may be provided 
by a burst of intense star formation, with implied star for- 
mation rates M* w 10 2 — 10 3 M Q yr~ 1 , or even higher for 
HyLIGs. At the high end of this interval, a starburst lasting a 
typical 10 8 yrs would produce a stellar mass of ~ 1O 11 M0, 
comparable to the luminous mass of a typical galaxy. There- 
fore, if powered only by star formation, the HyLIGs form stars 
at the rate expected for high redshift galaxies in their forma- 
tion process, undergoing an initial starburst that builds up the 
bulk of their stellar population. AGNs are known to be present 
in most of the known HyLIGs, but since the fraction of the 
FIR luminosity that they provide is unknown, the presence of 
intense star formation is not ruled out. Indeed, several lines 
of evidence indicate the presence of vigorous star formation 
in IRAS F10214 + 4724 (Rowan-Robinson et al. 1993; Kro- 
ker et al. 1996; Green & Rowan-Robinson 1996), in addition 
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to an embedded AGN (Elston et al. 1994; Soifer et al. 1995; 
Goodrich et al. 1996). Furthermore, since most of the known 
HyLIGs have been found in surveys of AGNs, the ubiquity of 
AGNs in known HyLIGs may be entirely a selection effect. 
Thus an unbiased survey for HyLIGs based on a sample of faint 
far-IR sources would be valuable for a reliable assessment of 
the nature of these objects. In addition, since HyLIGs can be 
observed to very significant redshifts, a determination of their 
number densities will strongly constrain the evolution of IRAS 
galaxies at these redshifts, as far as the most luminous part of 
the luminosity function is concerned. 

In this paper we present a survey for the most luminous 
galaxies detected by IRAS in a 1079 dcg 2 area of sky. We dis- 
cuss the sample selection process in Sect. 2 and our observation 
and reduction procedures in Sect. 3. The results are presented 
and discussed in Sects. 4 and 5, while our conclusions are sum- 
marized in Sect. 6. 

2. Sample selection 

We based our survey on a sample of sources selected from 
the IRAS Faint Source Catalog (FSC; Moshir et al. 1992). 
Since our project involves the optical identification of faint 
sources, it is important to avoid spurious FSC sources (result- 
ing from e.g., small-scale structure in 60 /im cirrus). We there- 
fore selected a survey area where diffuse 60 /im emission is 
faint (cf., the IRAS 60 /im maps presented by Rowan-Robinson 
et al. (1991b)). The area selected consists of the 4 hour 
R.A. (B 1950.0) interval between 21 h and l h , at Dec. (B 1950.0) 
less than —30° and Galactic latitude less than —40° We also 
used a 60 /im flux cutoff of 0.2 Jy, since below this value the 
FSC becomes rapidly less complete. The region selected con- 
tains 3057 60 /im FSC sources brighter than 0.2 Jy over an area 
of 1079 deg 2 . Stars and nearby galaxies were rejected by ex- 
cluding all sources detected at 12 /im, leaving 2719 objects in 
the sample. As a further safeguard against spurious sources, we 
used the FSC flux quality indicators and cirrus flag, to retain in 
the sample only those sources with a high-quality detection at 
60 /im and no confusion by cirrus. Since the spectral energy 
distribution (SED) of ULIGs peaks in the rest-frame 60 /im re- 
gion, ULIGs at z > 0.3 will have flux densities rising mono- 
tonically with wavelength in the IRAS bands. Therefore the 
resulting sample was further reduced by retaining only sources 
detected at both 60 and 100 /im. However, following Clements 
et al. (1996), sources with Sino/S^o > 5 (where Sioo and Sqo 
are the 100 and 60 /im flux densities as given in the FSC) were 
excluded, since such cold sources most likely arise from small- 
scale structure in Galactic cirrus. Finally, we rejected sources 
with associations in other catalogs as indicated in the FSC, 
thus limiting our sample to 313 objects. As shown in Sect. 5, 
these strict selection criteria make our survey a sparse (approx- 
imately 1 in 8) but unbiased survey for infrared galaxies with 
5*60 > 0.2 Jy over the 1079 deg 2 survey area. 

In order to select from our sample the most distant objects, 
we define the FIR loudness R by 

R^L m /L B = S 60 10 OA{B - liA5 \ (1) 



(see Clements et al. 1996), where Lgo is the 60 /im luminosity, 
Lb the luminosity in the B band, B the £>-band magnitude, and 
S 60 the 60 /im flux density in Jy. The bivariate B-60 /im lumi- 
nosity function has been derived by Saunders et al. (1990), who 
show that R increases monotonically with L e0 - This depen- 
dence accounts for the fact that the high luminosity cutoff of the 
luminosity function is much sharper in the optical regime than 
in the infrared. Therefore, 5 60 can be combined with the appar- 
ent B band magnitude to calculate R and hence obtain a crude 
estimate of the far-IR luminosity and distance of the object. 
An approximate B magnitude of the most luminous sources 
in our sample can be estimated as follows. For the cosmolog- 
ical parameters adopted here, a ULIG with Lir = 10 12 L Q 
will have Seo = 0.2 Jy if it is at z « 0.35. Using the bi- 
variate luminosity function of Saunders et al. (1990), 95% of 
these will have R > 10, and they will have a mean absolute 
B magnitude of — 20^0, or an apparent magnitude B w 21™0 
at z « 0.35. At these magnitudes, sources can be identified on 
optical Schmidt survey plates. Furthermore, since the IRAS er- 
ror ellipse for our sample sources is typically 10" x 30", and 
the extragalactic source density at B < 21 m is about 1000 per 
square degree, there is only about a 2% probability of chance 
superpositions at these magnitudes. Therefore, ULIGs in the 
FSC can be identified in optical Schmidt surveys and selected 
based on their Lqq/Lb ratio. This method has been succesfully 
used by Clements et al. (1996), to find 91 ULIGs with a me- 
dian redshift between 0.2 and 0.3, and a maximum redshift of 
0.43, by selecting only R > 10 objects from FSC sources iden- 
tified on Schmidt plates. This reasoning suggests that HyLIGs 
could be found in the same way, but selecting only sources with 
R > 100 (e.g., IRAS F10214+4724 has R = 350). However, 
a HyLIG with L m = 10 13 L Q will have S 60 = 0.2 Jy at z « 1, 
and a most likely B « 23™7. Such sources are below the plate 
limit of common Schmidt surveys, while with deeper imag- 
ing the density of faint sources becomes so high, that reliable 
identification is no longer possible without additional informa- 
tion. In order to circumvent these problems we have first ob- 
tained a sample of candidate distant HyLIGs by selecting those 
sources for which no reliable identification can be found on op- 
tical Schmidt plates, or which have very faint optical counter- 
parts. Optical follow-up (see Sect. 3) was then used to obtain 
the correct identifications and redshifts for the candidates. We 
note that the existence of faint but reliable FSC sources with- 
out optical counterparts above the typically Bj = 22 m limit 
of Schmidt surveys has been noted by several groups (Wols- 
tencroft et al. 1986; Rowan-Robinson 1991; Sutherland et al. 
1991;Clements et al. 1996; Oliver et al. 1996). Our programme 
is the first published project to systematically identify these op- 
tically faint and potentially very distant sources. 

We carried out an identification programme for our FSC 
sources on the U.K. Schmidt Telescope southern sky survey 
plates, digitized using the COSMOS plate scanning machine 
(Yentis et al. 1992). The COSMOS catalog provides Bj mag- 
nitudes to a completeness limit of Bj w 22 m , positions, major 
and minor axis lengths, position angles, and for objects with 
Bj < 21 m a classification as star or galaxy, based on an algo- 
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rithm described by Heydon-Dumbleton et al. (1989). The iden- 
tifications were performed using the likelihood ratio method 
(see Sutherland & Saunders (1992) for a detailed discussion). 
Briefly, the method assigns to every optical source a likelihood 
ratio 

e -r 2 /2 

£ _ (2) 

2-Kcna2N{<Bj), 

where r is the distance between the FSC and optical positions 
in a coordinate system where the IRAS error ellipse is a circle 
of unity radius, 




In these expressions, u\ and cr 2 are the major and minor axes 
of the FSC error ellipse, di and c?2 are the position differences 
of the optical source with respect to the FSC source projected 
on these axes, and N(<Bj) is the density of objects brighter 
than the candidate object. 

In calculating L, it is important to take into account 
possible errors in the star/galaxy classification performed by 
COSMOS. We noted that a number of FSC sources had coun- 
terparts with a very high value of L, which were however clas- 
sified by COSMOS as stellar, but with axial ratios significantly 
exceeding unity. During our observing programme described 
in Sect. 3 we obtained B-band images of 19 of these, covering 
a range of magnitudes and axial ratios. In these 19 fields, we 
found that all objects classified as stellar by COSMOS, were 
in fact galaxies if they had Bj < 18 m and axial ratio ex- 
ceeding 1.27. Some objects with lower axial ratios were also 
misclassified as stellar. We therefore reclassified all objects 
with Bj < 18 and axial ratio greater than 1.27 as galaxies. 
In calculating L for objects classified as galaxies, we com- 
puted N(<Bj) taking into account only objects having the 
same classification. For objects classified as stellar or having 
Bj > 20™5 (making them too faint for useful classification), 
the calculation of N(<Bj) took into account all sources, re- 
gardless of classification. This method allows for the possibil- 
ity of misclassification, while still somewhat favouring objects 
classified as galaxies. 

In Eq. (2) the simplifying assumption is made that the po- 
sition errors in the FSC are Gaussian. This assumption is ap- 
proximately correct for small r, but the FSC position error 
distribution has wings which are stronger than Gaussian ones 
(Sutherland & Saunders 1992; Clements et al. 1996; Bertin et 
al. 1997). In order to take these wings into account, all optical 
sources classified as galaxies were assigned L = 5 if they were 
within 1' from the FSC position and had Bj < 19 m . 

The identification process consisted of calculating L as de- 
scribed above for every optical object within 2' of the FSC po- 
sition. Following Clements et al. (1996), we consider an op- 
tical identification reliable for L > 5. Of our sample of 313 
FSC sources, 302 had identifications with L > 5. Of the re- 
maining 11, 5 were found to have Bj < 19™5 objects within 
1.'25 of the FSC positions, which present plausible identifica- 
tions given the non-Gaussian wings of the FSC position error 



distribution. The remaining 6 had no plausible optical counter- 
part with Bj < 20™5. The best "identifications" for these FSC 
sources had L < 2. These 6 sources thus form our sample of 
candidate HyLIGs. We note that IRAS F10214+4724, which 
has B = 22™5 and is located outside the IRAS FSC error el- 
lipse, would also have been selected by this method, if it was 
located within our survey area. 

3. Observations and reduction 

In order to obtain identifications and redshifts of our FSC 
sources we used the ESO Faint Object Spectrograph and Cam- 
era (EFOSC; Buzzoni et al. 1984) on the 3.6 m telescope of the 
European Southern Observatory at La Silla, Chile, on the nights 
of September 6-8, 1994. Conditions were clear but not photo- 
metric and the seeing (measured in Gunn i-band) was typically 
l'/5. The detector was a 512 2 anti-reflection coated, thinned, 
back-illuminated Tektronix CCD with a pixel size of 0'.'61 in 
imaging mode. Our observing strategy was as follows. First, a 
short (1 — 2min) exposure of the field was taken in imaging 
mode using a £?-band filter (and for some objects also in R 
and/or Gunn i), allowing the detection of objects several mag- 
nitudes below the COSMOS plate limit. Subsequently, long-slit 
spectra were taken of potential counterpart galaxies, in order of 
decreasing likelihood ratio L, and including galaxies above and 
below the COSMOS plate limit, until an emission line galaxy 
was found. The spectra were taken with integration times of 10 
to 30 min, using the B300 and R300 grisms with a l'/5 slit, pro- 
viding a spectral resolution A/ A A « 300 from 3640 to 6860 A 
(B300) or from 5970 to 9770 A (R300). Wavelength calibration 
was derived from exposures of a HeAr lamp. Photometric and 
spectroscopic calibration was achieved by observations of the 
spectrophotometric standard L870— 2. Flatfields were obtained 
from lamp exposures. Data reduction was performed using the 
standard long-slit reduction procedures as implemented in the 
IRAF package. 

4. Results 

Our results are summarized in Table 1 . Positions refer to the 
positions provided by the COSMOS catalog, or for sources be- 
low the COSMOS plate limit to detections on our EFOSC im- 
ages, adopting the astrometry of the COSMOS plates. Position 
errors are about 1" r.m.s. Magnitudes Bj are values provided 
by COSMOS, except for IRAS F21243-4501, where magni- 
tudes are derived from our EFOSC imaging, calibrated using 
the COSMOS Bj magnitudes of other sources in the field. 
Since this procedure ignores the differences between Bj and 
the Bessel B filter used in EFOSC, Bj for this object is only 
approximate. Far-IR luminosities Lfir have been calculated 
from Lfir = 47t£)£Ffir where Z?l is the luminosity distance 
and 

— 2 = 1.8 x 10- 11 ( 2.58^ + (4) 

erg s 1 cm 1 \ Jy Jy / 

(see Sanders & Mirabel 1996 and references therein). This 
procedure somewhat underestimes the total far-IR luminos- 
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Table 1. Parameters of distant FSC sample sources 



name R.A. Dec. z Seo >Sioo Bj £fir R 

(B 1950.0) [Jy] [Jy] [L ] 



IRASF00320- 


3307 


00 h 


32 m 


00?99 


-33° 07' 44'.'6 


0.439 


0.43 


0.87 


21 m 10 


4.9 x 10 12 


200 


IRASF00417- 


-3358 


00 h 


41 m 


40?41 


-33° 58' 03'.'0 


0.461 


0.24 


0.59 


21 m 97 


3.3 x 10 12 


240 


IRASF21065- 


3451 


21 h 


06 m 


33?93 


-34° 51' 34'.'0 


0.329 


0.36 


1.55 


21 m 54 


3.1 x 10 12 


250 


IRASF21243- 


4501 


21 h 


24 m 


26?05 


-45° 01' 54'.'2 


0.834 


0.30 


0.90 


23 m 7 


1.9 x 10 13 


1500 


IRASF22148- 


4013 


22 h 


14 m 


52f07 


-40° 12' 58'.'4 


0.529 


0.33 


0.83 


21 m 6 


4.4 x 10 12 


240 


or: 




22 h 


14 m 


51?96 


-40° 13' 02'.'7 


0.380 






21 m 5 


2.0 x 10 12 


220 


IRASF23555- 


3436 


23 h 


55 m 


32?14 


-34° 36' 29'.'3 


0.490 


0.31 


0.71 


20?98 


4.8 x 10 12 


130 



ity because it does not include a if -correction. An accurate 
if -correction is not possible because of the lack of knowledge 
of the SED of the sources. However, under the assumption that 
the SED is similar to that of the prototypical ULIG Arp 220, 
we find that the underestimate introduced by Eq. (4) could be 
up to 50% for the most distant objects. 
Notes on individual sample sources: 

- IRAS F00320-3307: while classified as only one galaxy 
by COSMOS, this system consists of 2 interacting galaxies 
at z = 0.439. The compound spectrum shows strong [O II] 
3727 A, in addition to [Nem] 3869 A, [Om] 5007 A and 
Ha. 

- IRAS F00417-3358: the object with the highest likelihood 
ratio in this field, and therefore the a priori most likely 
counterpart, was found to be a luminous object showing 
[O II], [Ne III], Ca H and K absorption, and a 4000 A break 
at z = 0.461. 

- IRAS F2 1065 -3451: the second most likely counterpart 
as indicated by our identification process, barely resolved 
at z = 0.329, and showing strong [Oil], in addition to H/3 
and [Om]. 

- IRAS F21243-4501: none of the possible COSMOS iden- 
tifications showed emission lines, but a fainter object close 
to the IRAS error ellipse was found to have [Ne v] 3426 A 
and [Oil] at z = 0.834. A broad feature at the expected 
wavelength of Mg II 2798 A may also be present. 

- IRASF22148-4013: spectra of two galaxies only 4'/5 
apart yielded redshifts of 0.380 and 0.529, based on strong 
[O II], H/3 and [O III] lines (for both objects) and also strong 
Ha, [Nil] 6584 A and [S II] 6716 and 6731 A lines (in the 
object at z = 0.380). The two objects are of closely similar 
Bj magnitude. Either of these may be the correct identi- 
cation, or they may both contribute part of the FSC 60 yum 
flux density. In either case at least one of the objects is a 
ULIG, but none is a HyLIG. 

- IRAS F23555-3436: the most likely identication from the 
COSMOS plate is a distorted object outside but close to the 
FSC error ellipse showing [Oil] and [Nem] at z = 0.490. 

In addition, we observed one object not in our sample of 
6 candidate distant objects, IRAS F22569-5523, in order to 



check possible misidentification, since the only likely coun- 
terpart was classified by COSMOS as a fairly bright star with 
low axial ratio. However, this object is in fact a Bj = 17™86 
galaxy with a possible tail or extension towards the east, and 
showing [Nem], [Oil], H/3 and [Om] emission lines at z = 
0.235. The COSMOS position for this object is R.A. (1950) = 
22 h 56 m 53?73, Dec. (B1950) = -55° 23' 24'.'6 and its far-IR 
luminosity is 8.8 x 10 11 L Q . Accounting for the flux beyond 
100 ,um, this object is also a ULIG. 



5. Discussion 

All objects from our sample are found to have a FIR loudness 
R > 100. In contrast, the highest value of R among the L > 5 
sources that were removed from the sample is 76. Therefore 
our approach of selecting those sources which do not have reli- 
able counterparts above the COSMOS plate limit, or for which 
the counterpart is so faint that misidentification is no longer 
unlikely, proves to be very effective in selecting sources with 
extreme values of R. What is the nature of these objects? Of 
our sample of 6 sources, one is a HyLIG and five are non- 
hyperluminous ULIGs. The five non-hyperluminous ULIGs are 
all detected on the COSMOS plates and have Bj < 22.0 and 
R < 250. Their mean redshift z = 0.45 is higher than the 
highest known redshift of any non-hyperluminous ULIG prior 
to this study, indicating that our procedure is also a powerful 
method for selecting distant ULIGs. The HyLIG in our sample 
is the only object not detected on the COSMOS plates and this 
object has B = 23™7 and R = 1500. This result confirms that 
HyLIGs can be found by selecting objects with extreme values 
of R. The main difficulty in applying this method is the large 
size of the IRAS position error ellipses, which precludes a di- 
rect optical identification at the faint magnitude levels expected 
for distant HyLIGs. However, future surveys, such as the on- 
going European Large Area Infrared Survey (ELAIS; Oliver 
1996), and surveys with SIRTF and FIRST, and with SCUBA 
on the James Clerk Maxwell Telescope (JCMT) will provide 
substantially better positional accuracy and not suffer from this 
identification ambiguity. The method used here for selecting 
the most luminous and distant objects can be adapted directly 
to those surveys. 
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The small size of our sample, which contains only one 
HyLIG, precludes any detailed statistical inferences, which 
must await more extensive programmes using this selection 
and identification method, based on IRAS data or on the sur- 
veys mentioned previously. However, a number of trends in 
our data merit further discussion. In the first place, the de- 
tection of [Ne v] emission in the only HyLIG in our sample 
shows that this object contains an AGN. Thus all three IRAS- 
selected HyLIGs discovered so far (IRAS F10214+4724, 
IRASF15307+3252 (Cutri et al. 1994; Hines et al. 1995) and 
IRASF21243-4501 (this work)) contain AGNs. While the 
statistics for HyLIGs is still based on small numbers, the re- 
sult is significant, since the [Ne v] line was not detected in 
any of the non-hyperluminous ULIGs in our sample, while our 
spectra did cover the wavelength where this line would be ex- 
pected. Thus the HyLIGs form a remarkable contrast with the 
non-hyperluminous ULIGs, where the presence or absence of 
AGNs is a strongly debated issue, and direct evidence for the 
presence for an AGN is very scarce. 

Our procedure brings about incompleteness in our sample 
of R > 100 objects in two ways: identification incomplete- 
ness and selection incompleteness. The former effect arises if 
objects with R > 100 fail to be selected by our L < 5 crite- 
rion, which occurs if a bright galaxy lies close to the line-of- 
sight to a distant FSC source, giving rise to erroneous identi- 
fication with the bright galaxy. As noted in Sect. 2, the prob- 
ability of misidentification in this situation is only about 2% 
for galaxies with Bj < 21 m . Since the large majority of our 
L > 5 identifications have counterparts significantly brighter 
than Bj = 21 m (for 85% of the objects with L > 5, the coun- 
terpart has Bj < 19™0), the probability of chance superposi- 
tions is much less than 2%, and the identification incomplete- 
ness can thus be neglected. 

However, the sample of 313 objects used for our iden- 
tification programme does suffer from selection incomplete- 
ness. Our selection method was aimed at rejecting spurious 
sources; however, as shown below, it must have removed a sig- 
nificant number of real sources from the sample as well. The 
relevant selection criteria are the requirement to have a high- 
quality 60 /im detection, no cirrus confusion, and a detection at 
100 /mi. While these criteria were effective at rejecting spuri- 
ous detections, they also introduce a selection incompleteness, 
and may have rejected some distant objects. In order to assess 
the magnitude of this effect, we compare our sample to the 
FSS-z I sample described by Oliver et al. (1996). This sam- 
ple has been constructed using low-cirrus regions with good 
IRAS 60 /mi coverage and is estimated to be 99% complete for 
Sao > 0.2 Jy, which is the same flux limit as the sample de- 
scribed in the present paper. It contains 1931 IRAS FSC galax- 
ies over an area of 839 deg 2 , giving a source density of 2.30 
per deg 2 . Adopting this source density as characteristic for the 
present survey shows that a total of 2483 expected IRAS FSC 
galaxies over the entire survey area should be expected, a plau- 
sible number given that, including spurious sources, our ini- 
tial extragalactic sample in this area contained 2719 objects 
(see Sect. 2). In contrast, only 313 objects were retained in 



our sample of candidate objects after the strict selection cri- 
teria described in Sect. 2 had been applied. However, since 
none of these criteria introduces a bias in luminosity or dis- 
tance, our sample is unbiased and our survey thus constitutes a 
sparse (approximately 1 in 8) survey of infrared galaxies with 
Seo > 0.2 Jy over the 1079 deg 2 area. Hence we can use our 
results to estimate a number density for HyLIGs at z < 1 of ap- 
proximately 7 x 10~ 3 deg -2 , with considerable uncertainty due 
to the small numbers involved. We note that, adopting the local 
60 /im luminosity function of Saunders et al. (1990), this esti- 
mate implies significant evolution in the infrared galaxy pop- 
ulation to z = 1. Only in the unlikely case that the HyLIG 
detected in our sparse survey was the only z < 1 HyLIG in the 
entire 1079 deg 2 survey area, no evolution would be needed. 

We finally note that since we are using L 60 /Lb to select 
luminous objects, our selection method is robust against the 
presence of gravitational lensing, provided the corresponding 
magnification factors are similar at 60 /mi and B. As a result, 
once a redshift and hence an infrared luminosity is available, R 
and Z/tr, may be combined to address the possibility of gravita- 
tional lensing. We illustrate the method using the lensed HyLIG 
IRASF10214+4724 and the HyLIG IRAS F2 1243 -4501, 
identified in the present work. As noted in Sect. 2, 
IRASF10214+4724 has R = 350. Using the bivariate B- 
60 /mi luminosity function of Saunders et al. (1990), we find 
a most likely intrinsic infrared luminosity L^ r of about 3 x 
10 12 L Q . The apparent luminosity following from the redshift 
of 2.28 on the other hand, is L^ p = 2 x 10 14 L . The large 
discrepancy between L\^ T and suggests gravitational am- 
plification by a factor of about 60. Using the same reasoning, 
for IRASF21243-4501 we find L£f = 1.6 x 1O 13 L and 
^ir P = 1-9 x 10 13 L©. Because of the similarity of the two val- 
ues, there is in this case no indication for gravitational lensing. 
Caution is required when applying this method, since the un- 
derlying assumption of similar magnification factors at optical 
and infrared wavelengths may easily be violated, as is the case 
in IRAS F10214+4724, where an optical magnification by ap- 
proximately a factor of 100 is found (Eisenhardt et al. 1996), 
whereas the infrared magnification is only approximately a fac- 
tor of 10 (Downes et al. 1995; Green & Rowan-Robinson 1996; 
Serjeant et al. 1998). Therefore the actual presence or absence 
of gravitational amplification must always be established by ad- 
ditional observations. However this method may be useful for 
selecting candidate gravitationally lensed sources for further 
study. 

6. Conclusions 

1 . We have identified the most luminous infrared galaxies in 
an unbiased sample of 313 reliable extragalactic IRAS FSC 
sources with S 60 > 0.2 Jy. Our method is based on the 
bivariate B-60 /im luminosity function of infrared galax- 
ies, which implies that the most luminous objects have the 
highest values of R (as defined by Eq. (1)), and have op- 
tical counterparts that are so faint that they cannot be reli- 
ably identified (or are undetected) in typical Schmidt sur- 
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veys. Using optical spectroscopy, we have systematically 
identified the optical counterparts of all of the 6 sources 
in our 60 fim sample that were too faint in Bj to be re- 
liably identified on the UKST plates. Our results confirm 
that this method selects the galaxies with the largest values 
of R, so that these galaxies are indeed the 6 most lumi- 
nous infrared galaxies in our sample. Five of these are non- 
hyperluminous ULIGs with a mean redshift of 0.45, higher 
than any previously known non-hyperluminous ULIG; the 
remaining source is a HyLIG at z = 0.834. 

2. The HyLIG in our sample (IRAS F21243-4501) contains 
an AGN, as shown by the presence of [Nev] emission. 
Hence all infrared-selected HyLIGs discovered so far un- 
ambiguously show the presence of AGNs. In contrast, none 
of the non-hyperluminous ULIGs in our sample show evi- 
dence for the presence of AGNs, and such evidence is rare 
among non-hyperluminous ULIGs in general. 

3. Our method is robust against the effects of gravitational 
lensing if the optical and infrared magnification factors are 
similar. Under this assumption this method may be use- 
ful for selecting candidate gravitationally lensed sources by 
comparing an intrinsic luminosity (estimated from R) with 
the apparent luminosity (calculated from S^o and z). 

4. Our survey consitutes an unbiased, sparse (approximately 
1 in 8) survey of infrared galaxies with Sea > 0.2 Jy over 
a 1079 deg 2 area, and the results allow an estimate of the 
number density of HyLIGs at z < 1 of approximately 
7 x 10~ 3 deg~ 2 , with considerable uncertainty due to the 
small numbers involved. Compared to the local luminosity 
function of infrared galaxies, this estimate indicates sub- 
stantial evolution at the highest luminosities, except in the 
unlikely case that the HyLIG found in our sparse survey is 
the only HyLIG at z < 1 in the entire 1079 deg 2 survey 
area. 
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